function firTest()
    % Load the reference data
    data = load("../data/gse2.txt");
    % Set an FIR filter: fir1(50, 0.4); %Uses Hamming window
    b = [-0.000000000000000, -0.001056235801065,  -0.000769341020163, ...
          0.000956323223723,  0.001976082742122,  -0.000000000000000, ...
         -0.003265384800345, -0.002568519852901,   0.003234633130890, ...
          0.006519908075213, -0.000000000000000,  -0.009825739114867, ...
         -0.007365685405410,  0.008881348924986,   0.017256056989442, ...
         -0.000000000000000, -0.024784271698734,  -0.018417666768131, ...
          0.022299534288278,  0.044222443880910,  -0.000000000000000, ...
         -0.071469809226860, -0.060430328816090,   0.092317626953209, ...
          0.302027315266443,  0.400523418058701,   0.302027315266443, ...
          0.092317626953209, -0.060430328816090,  -0.071469809226860, ...
         -0.000000000000000,  0.044222443880910,   0.022299534288278, ...
         -0.018417666768131, -0.024784271698734,  -0.000000000000000, ...
          0.017256056989442,  0.008881348924986,  -0.007365685405410, ...
         -0.009825739114867, -0.000000000000000,   0.006519908075213, ...
          0.003234633130890, -0.002568519852901,  -0.003265384800345, ...
         -0.000000000000000,  0.001976082742122,   0.000956323223723, ...
         -0.000769341020163, -0.001056235801065,  -0.000000000000000];
    % Filter data
    yfilt = filter(b, [1.0], data);
    % Write the data
    fileID = fopen('firReference.txt','w');
    for i=1:length(yfilt)
      fprintf(fileID, '%.13e\n', yfilt(i));
    end;
    fclose(fileID);
end
